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Abstract 


«D 


Numerical and experimental simulation of un- 
steady airflow through the control vaJve and slot- 
ted air due" of a circulation control rotor is de- 
scribed. The numerical analysis involves the solu- 
tion of the quasi-one-dimensional compressible 
fluid-dynamic equations in the blade air duct to- 
gether with the coupled isentrop 4 c ti^v equations 
for flow into the blade through the valve and out 
of the blade through the Coanda slot. Numerical 
solutions are compared with basic experimental re- 
s’ Its obtained for a mockup 01 a circulation con- 
trol .otor and its pneumatic valving system. The 
pneumodynamic phenomena that were observed are dis- 
cussed with particular emphasis on the characteris- 
tic system time lags associated with the rerponse 
of the flow variablps to i oansient and periodic con- 
trol valve inputs. 


Notation 


A Cross-sectional area of duct 

Effeccive expansion area at valve exit 
A^ Valve area 

^ Duct friction factor 

C g Slot discharge coefficii it 

C v Valve discharge coefficient 

c v Specific heat at constant volume 

Up Hydraulic diameter of duct 
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e To cal energy per unit volume. v^ 

F Flux times area vector 
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g 

0 Vector of nonhomogeneous terms in flow 

equations x 


h Heat transfer coefficient 
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PRECEDING FACS BEANE NOT FttMEO 


Average duct height 
Thermal conductivity 
Length of duct 

Mass flow rate 

Mass flow rate per init area in duct 

Mass flow rate ’-er ui it area computed from 

isentropic flov ; neor/ 

Mass flow rate per unit urea through slot 

Mass flow rate per unit arec through valve 

Static pressure in duct 

External pressure at slot exit 

Total pressure i* ct 

Duct static tempt .are 

P ' • to ’ temperature 

Du 't wail temperature 

Time 

Vector of dependei w variables 

Vector of diffusion terms 

Average duct velocity in spanwise or x 
direction 

Approximate average duct velocity in cncr^- 
wiue or y direction 

Slot width 

Spanwise coordinate 

Coordinate at duct eutrance 





au£L*x.\ 
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i. 



x end Coordinate at duct end 

y Chordwise coordinate 

V Ratio of specific heats 

p Density 

w Angular frequency 

Subscripts 

I Variable computed from isentropic flow 

at valve exit 

pi Variable in pressure supply plenum or 

rotor hub 

Prefix 

Peak-to-peak value of variable 


The present paper reviews the modeling techniques 
described in Reference 1 and presents results ob- 
tained by applying the code to a simple experimen- 
tal mockup of a model rotor blade and its pneumatic 
valving system. The discussion of the results, 
from both the HFPA predictions and the experiment, 
focuses on the dynamic response of the flow vari- 
ables to periodic and transient control valve in- 
puts. These response characteristics, particular- 
ly with reference to the characteristic system 
time lags, have implications for the performance 
of actual rotor systems. 

Theory 

The modeling techni *es employed *n the pre- 
sent research are described in Reference 1 where 
emphasis is on the details of the theoretical 
formulation and numerical technique. The theory 
is presented here in abbreviated form; a more com- 
plete discussion is contained in Reference 1. 


T ntroduction 


Basic Equations 


Circulation control rotor (CCj technology 
applied to rotary-wing aircraft or stopped-rotor 
vertical takeoff and landing (VTOL) aircraft, such 
as the X-Wing, offers several distinct advantages 
over conventional rotor technology. For rotary- 
v '.ng aircraft, CC . echnology provides a solution 
to such problems as high vibration levels, retreat- 
ing blac*'' stall, hub/pylon drag, and implementing 
higuer harmonic control. For stopped-rotor air- 
-aft, the historical limitations of aeroelastic 
divergence, flutter, blade dynamic instability, 
and critical resonance conditions during rpm re- 
duction are eliminated. Both types of aircraft 
use a shaft-driven rotor with blades having circu- 
lation control airfoils which generate lift through 
the Coanda principle. The CC airfoils employ a 
rounded trailing edge with a thin jet of air tan- 
gentially ejected from a spanwise slot adjacent to 
the rounded (Coand^) surface. The jet of air sup- 
presses boundary 1 ver separation and moves the 
rear stagnation streamline toward the lower sur- 
face, thereby increasing lift. Lift is increased 
in proportion to the mass flow rate of compressed 
air in the jet. Pitch and roll control require- 
ments are obtained by cyclic modulati #n of the 
mass flow rate with valves in the nonrotating 
system. Higher harmonic cyclic control can be 
similarly applied for reducing blade stresses, 
transmitted shears, vibration, and power require- 
ments. 

The design of CC rotors and their pneumatic 
control systems requires an understanding and 
appreciation of the phenomenology involved in the 
control and distribution of airflow to the Coanda 
slots in these rotors. A capability for analyti- 
cal prediction of rotor Coanda airflow is also 
essential. The term "pneumodynamics" has been 
coined to refer to the aerodynamic response char- 
acteristics of CC rotor system internal airflow. 

In a recent paper, Watkins et al.^ describe 
the modeling techniques employed in the HFPA (high 
frequency pneumodynamic analysis) computer code 
developed for CC rotor pnetimod>namic analysis. 


The HFPA code solves the quasi-on* dimension- 
al, compressible fluid dynamic equations for un- 
steady flow in a spanwise blowing air supply duct 
internal to a rotating CC rotor blade, together 
with the coupled isentropic flow equations for 
flow into the blade through the control valve and 
out of the blade through the Coanda slot. The 
duct flow equations have che conservative form 
representation 

U, + t F, + G + V » 0 (1) 

* t A *x 

where the vector of dependent variables is 
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and the flux is given by 
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The vector of nonhomogeneous terms, which includes 
the effects blade rotation at angular velocity 
u, wall skin friction fact— i t (x), wall heat 
transfer coefficient h(x), and mass flux through 
the slot m s (x,t) is 
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For completeness, the vector of axial diffusion 
ter™s i*? included, althouph it is insignificant in 
the present problem. It is 

r o i 



Figure 1 shows the idealized representation adopted 
for the cam type valve of interest in the present 
study. The controlled value of the valve opening 
is designated as A^t), and A e is the fixed effec- 
tive opening at the duct entrance. Like the slot 
flow, one-dimensional isentropic flow theory with 
a discharge coefficient is used to represent flow 
through the valve. The valve mass flux is 



The equation of state is 


P 


2 




“(t) " c v m I <P 0 * T 0 »P (x e> t)) < 5 > 

Pi Pi 

where pQ and Tq are the total p\essure and tempera- 
ture, respectively, in the air supply plenum in the 
rotor hub and p(x e ,c) is the back pressure in the 
duct entrance downstream of the valve. When the 
(2) flow is choked, the dependency on p(x e ,t) is eli- 
n:5nated. 


In the above relations, the conservative variables 
ra and e are defined as 


The other boundary condition, at x =i x ent |, is 
the stagnation of the flow at the end of the duct. 


x 

2 2 
v v 

e « P (c v T + ) 

and w (x,t) and m s (x,t) are the local (with respect 
to spinwise x location) slot width and blowing mass 
flux, respectively. 


'Tie discharge coefficients C y and C s , used cs 
correlation parameters for a given system relate 
the geometric area of a system to the effective 
arec. For simple configurations such as orifice 
plates or venturi tubes in a well conditioned sys- 
tem where the airflow is well behaved, the discharge 
coefficient value can be obtained from a standard 
engineering handbook. 


The small average chordwise velocity component 
inside the duct can be crudely approximated as 


v y - 2 PHd 


(3) 


in which ^(x) is the local effective height of 
the duct. 


Hie mass flux distribution through the Coanda 
slot is assumed to be Wtll represented by one- 
dimensional isentropic flow theory with a discharge 
coefficient C g (x); 

V*» e > “ C s m I fp 0* T 0' P ext ) (4) 


The above functional relationship is the usual ex- 
pression for isentropic flow expanding through a 
nozzle from a plenum (in this case the duct in- 
terior) to ambient conditions (in this case the 
exterior of the rotor blade). p^(x,t) and T^(x,t) 
are the local (internal) total pressure and tem- 
perature, respectively , and Pe x t^ x,t ^ is external 
pressure a the exit of the Coanda slot. Depen- 
dency upon external pressure is eliminated when 
the flow is choked. 


Boundary Conditions 

At the upstream boundary x - x fi , a plenum 
Supplies blowLng air to the duct through a valving 
system. The valve opening area has constant and 
rotor-azimut.i-dependent components regulated by the 
collective and cyclic components, respectively, ^f 
the aircraft control system. The mass flow and 
duct pressures are thus determined by these valve 
settings in concert with the duct/slot configuration. 


Finite Difference Method 

It. applying the theory of the previous section 
to the analysis of a rotor blade, the entire blade, 
including the transition duct (ducting from hub 
valve exit to airfoil portion of rotor blade), is 
discretized by dividing it into spanwise (radial; 
segments. The governing differential equations, 

Eq. (*), are solved by the implicit, "delta form," 
finite-difference procedure of Beam and Warming^ at 
the grid points located at segment boundaries, in- 
cluding the valve exit as shown in Fig. 1 and the 
duct end. The solution of the duct flow is active- 
ly coupled \o the flow through the valve and through 
the upstream boundary condition. With the exception 
of the upstream boundary condition for flow into the 
duct from the valve, the numerical analysis is rather 
standard and is described fully in Reference 1. 


In Reference 1, expressions are presented for 
an approximate numerical upstream boundary condi- 
tion based on deriving for various valve types an 
auxiliary relationship by ignoring the time depen- 
dent term in the momentum equation between the 
first two grid points (grid point 1 and 2 in Fig. 
1), The approximation for the idealized valve of 
Fig. 1 is 


p l A e + 


(P 1 +P 2 ) 

~2 


(A 2”V “ p 2 A 2 + p 2 A 2 


(A e ~A 2 )p 2 u 


2 2 ( 'f m 2 

(Xj-X 2 ) - 2 ( Xl -x 2 ) 


( 6 ) 


Approximations such as in Eq. (6) were applied be- 
cause the investigators were unsuccessful in aply- 
ing the more standard, and less approximate, pro- 
cedures for inlet boundary conditi;.is of applying 
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one-sided differencing or the method of character- 
istics. The presence of the valve poses special 
difficulties in the application of these standard 
techniques. In a recent communication, Anderson of 
United Technologies Research Center imiieaLed that 
he was successful in applying the method of char- 
acteristics to derive inlet boundary conditions 
for the HFPA code for a sudden expansion or gate- 
type valve. The extension of his procedure to the 
cam valve type shown in Fig. 1, is not straight- 
forward due to the 90 deg change in flow direction. 
Since this valve is he type considered in the pre- 
sent investigation, Eq. (6) was retained. 

Numerical Solution proceuure 

Hie basic approach followed in obtaining nu- 
merical solutions to control inputs, consisting of 
constant (collective) settings or constant collec- 
tive settings with superimposed constant amplitude 
cyclic inputs, is to obtain the steady state or 
time periodic solutions by allowing them to evolve 
from a transient. 

The solution starts with the system essential- 
ly at supply pressure conditions. The initial 
valve position is assumed to be full/ open as speci- 
fied by the sura of its Fourier components supplied 
3S input. The solution is then advanced with time 
as the flow in the blade adjusts itself to reach 
steady state flow conditions for the fully open 
valve. If cyclic components are to be considered, 
the cyclic input is imposed after fully onen steady 
state flow conditions are achieved. The solution 
is thereafter advanced further with time until 
periodic conditions are attained. A harmonic 
analysis of pressures, temperatures, and mass flow 
rates is performed on the results from the final 
cycle. 

Experiment 

The analytical formulation and the numerical 
methods described in the previous sections were 
validated by applying them to a basic experiment 
and comparing the results with the experimental 
data^. The experiment was performed using a mock- 
up of a rotor b’ade pneumodynamic system incor- 
porating a Vc.^ve similar to the idealized valve of 
Fig. 1, The model consisted of a supply plenum 
and two rotating cams with different profiles 
mounted on the same shaft allowing either a one- 
per-rev or two-per-rev profile to be selected by 
axially positioning the cams on the shaft. The 
flow was modulated by means of the cams from the 
plenum through a single nozzle into a stationary 
plugged pipe with a slot. 

The pipe with a full-length slot represented 
a 3ubscale circulation control rotor blade. The 
experiment was originally performed^ to establish 
the feasibility of modulating the weight flow in a 
one-per-rev and two-per-rev manner and to show that 
the three components of the airflow were additive. 
For simplicity, the correct relative motion was es- 
tablished by interchanging the rotating ana non- 
rotating components. The plenum and pipe (hub and 
rotor blade) were stationary, and the cams (the 



Fig. 1 Idealized valve 

valving system) rotated. The pipe incorporated 
several pressure taps distributed along its length. 
The model is shown in Fig. 2, which is reproduced 
from Reference 1. 
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Fig, 2 Experimental model 


Dynamic measurements of pipe pressure were ob- 
tained for a range of cam rotational speeds. Ir. 
addition, quasi-dynamic results were obtaiied by 
positioning the cam at discrete increments of 30 
deg. Average mass flow data were collected from a 
venturi tube located upstream of the plenum. 

Figure 3, also from Reference 1, shows the 
measured valve area formed by the one-per-rev and 
»:wo-per-rev cams spaced at a minimum distance of 
0.01 in from the nozzle inlet. The points shown 
were calculated from the measured gap between the 
cam profile and the nozzle; the area variations are 
intended to approximate sinusoids as shown. The 
two-per-rev cam profile is a rather crude approxi- 
mation and contains a significant higher harmonic 
content. The measured values Wv rc used in the nu- 
merical predictions. 
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Results and Discussion 


Numerical results were obtained for the expert- 
rental configuration of Fig. 2 by discretizing the 
pipe into 16 segments of approximately equal length. 
The pipe was assumed to be adiabatic, and the fric- 
tion coefficient was represented by the formulas 
for fully developeu flow. The integration time- 
step for the calculations simulating the quasi- 
dynamic experiments was 5x10”^ ^fec. For the dy- 
namic calculations, valve cycling was imposed 
after an elapsed time of 0,25 cec. and the time- 
step was reduced to 1200 steps per revolution there- 
after. In general, the results were qualitatively 
similar for the one- and two-per-rev cams. There- 
fore, the results presented here are, with one ex- 
ception, for the one-per-rev cam. 

Typical quasi-dynamic numerical and experimen- 
tal results, indicating the dependency of mass flow 
rates on nozzle control area at a constant slot 
height of 0.042 in., are displayed for the one-per- 
rev cam in Fig. 4. In this figure, the mass flow 
rate is presented as a function of total valve area 
at a set plenum pressure, *.a the valve area in- 
creases, the data show that mass flow rate tends to- 
ward a constant value, indicating that the flow has 
become controlled by the slot opening. 



Fig. 4 Mass flow variation 
Mass Flow and Pressure Comparison 

Figure 4 shows the sensitivity of the discharge 
coefficients C v and C t for a wide range of valve 
areas. Also shown is the effect of assuming con- 
stant values of 0.95 and 0./5 for both the valve and 
slot discharge coefficients in the numerical solu- 
tion. The results obtained for the higher value 
Cy * C s ■ 0.95 agree more favorably with the experi- 
mental results, although the lower discharge coe : f i- 
cient is likely to be realistic physically. That a 
concept as simple as a discharge coefficient can be 
used to correlate the mass flow and pressure losses 
is very encouraging; this is discussed in ^eoth in 
Reference 4. For the remainder of the theoretical 
results 0^ « C s ** 0,^5 were used without trying for 
a better correlation. As indicated in Reference 1, 
strict interpretation of Cy and C s is probably not 
desirable since the various deficiencies and approxi- 
mations in the computer model can, to some degree, 
be absorbed by adjusting them. Moreover, the con- 
cept of a constant nozzle discharge coefficient is 
ualy a first approximation for the actual flow loss- 
es which occur over a range of mass flow rates in 
the analysis. Recent experience with the HFPA code 
indicates that more physically realistic discharge 
coefficients can be obtained at the expense of de- 
creased resolution of wave reflection phenomena by 
including second-order numerical damping in the 
HFPA algorithm. 

Figure 5 illustrates the numerical and experi- 
mental cyclic variation of total pressure for a 
typical dynamic case. The numerical results of 
T ig. 5 were obtained with * C g - 0.95, The nu- 
merically computed pressure profiles have an 
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azimuthal phase lag from the valve area profiles of 
Fig. 3 of an approximately 35 deg phase shift. For 
this configuration, the phase shift is principally 
due to the finite speed of wave propogation, and 
is referred to as "sonic lag." A smaller portion 
of the delay is due to the "capacitance lag" effect 
caused b; the finite pipe volume. The phase dif- 
ference between the valve area setting and the pres- 
sure response to it has been eliminated from Fig. 

5 because the experimental phase angle results were 
not accurate. 



Fig. 5 Total pressure variation 

Table 1 contains data obtained from a harmonic 
analysis of the dynamic numerical results of the 
case for which the pressures are plotted in Fig. 5, 
The numerical phase lag has been retained. A har- 
monic analyclo of the experimental pressure data 
obtained from the experimental curve presented in 
Fig. 5 is also shown for comparison. Magnitude and 
phase for the valve opening area, total mass flow 
rate, and total pressures at two locations are given 
from tne numerical results where the phases are re- 
ferred to the maximum valve opening at zero deg. 


(The experimental pressure phase shifts given in 
the table have an arbitrary reference). The magni- 
tudes are all normalized with respect to the peak- 
to-peak values. 

Table 1 shows that the summation of the higher 
harmonic content of the mass flow rate output is a 
rather significant 20 percent while the one-pec-rev 
valve control input contains only approximately 10 
percent. Therefore, the mass flow exhibits some 
nonlinear amplification. The higher harmonic con- 
teat of the pressure response does not exhibit si- 
milar amplification. 

The phase shift phenomenon is clearly illustra- 
ted by the data in the table. The first harmonic 
of the mass flow rate lags the valve opening by 39 
deg. The first harmonic pressure at, x/L = 0.524, 
lags the valve opening by 3t> deg. At x/L - 0.804, 
this lag increases to 42 deg. due to the finite 
wave propogation speed or sonic lag. Because of 
:he specification of a zero numerical integration 
dissipation parameter in the numerical algorithm 
used for the eye ic calculations, computed results 
cental:: slight, stable spatial numerical oscilla- 
tions which make it difficult to place confidence 
in the exact amount of the lag in pressure predic- 
ted between locations. Baned on the speed cf sound, 
the pure sonic lag between the two points shown 
should amount to only about one-half the lag shown. 
The remali ing lag uay be due to a physical phenomenon 
or to numer 5 cal uncertainty. Although, as discussed 
previously, tu® problem of spatial oscillations can 
be reduced by employing numerical damping, this 
damping suppresses rr^tiple wave reflections and 
renders the code less aole to predict resonance 
phenomena. Fortunately. :he computed mass flow 
rate is the result of a spatial integration which 
has the effecc of smoothing the spatial oscillations 
and should be fairly accurate. 

Peak-to-Feak Pressure Comparison 

Figure 6 summarizes the results from the dyna- 
mic calculations in the form of peak-to-peak total 
pressures versus frequency for the one-per-rev can. 


Table 1, Harmonic analysis of typical results 


a 

|. 

A n /dA 

(dA-0.8904 IN 2 ! 
MAG {PHASE, DEG) 

m n /dm 

(dm-0, 00331 SlUG/SEC) 
MAG (PHASE, DEG) 

POn^PO 
At x/L * 0.804 

°0*. /d P0 
At x/L « 0.524 

COMPUTATION 
dpo-3.71 LB/IN 2 

MAG(PHASE.DEG) 

EXPERIMENT 
dpg*3.99 LB/IN 2 
(APPROXIMATE) 
MAG(PHASE.DEG) 

COMPUTATION 
dpg-3.77 LB/IN 2 

MAGtPHASE.DEG) 

EXPERIMENT 
dpo-3.76 LB/IN 2 
(APPROXIMATE) 
MAG(PHASE.DEG) 

0 

0.514 

0.686 

0.425 

0.619 

0.432 

0.637 

1 

0.509 (353) 

0.488 (32) 

0.520 :35) 

0.404 (360) 

0.523 (29) 

0.500 (350) 

2 

0.C27 (197) 

0.086 (241) 

0.016 (55) 

0.034 (176) 

0.017 (131) 

0.043 (171) 

3 

0.016 (128) 

1 

0.0^4 (100) 

0.031 (262) 

0.019 (205) 

0.007 (276) 

0.021 (203) 

4 

1 

0.006 (26! 

0.005 (2<;| 

0 002 (179) 

0.007 15) 

0.003 (78) 

0.005 (45) 

6 

0.004 (322) | 

0.002 (63) 

0.002 (179) 

0.01 i (347) 

0.003 (346) 

0.007 (30) 

| MOTE: w/frt * 49.7 Hj AND p p f - 8.89 LB/IN 3 . 
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In Fig. 6 the numerical solution is able to pre- 
dict fairly well the trends and magnitudes of the 
pressure response. One possible source of error 
in the numerical calculations is the assumption of 
constant discharge coefficient for the wide range 
of area variation. Somewhat more remote are the 
possibilities for experimental error. 



Fig. 6 Peak-to-peak total pressure 

Figure 7 plots the dynamic results of the peak- 
tu-peak total pressure near the beginning of the 
slotted portion of the pipe against plenum pressure 
for the one-per-rev cam at two different frequencies. 
The agreement between the measured and numerical re- 
sults in this figure, in general, is seen to be 
fair. As shown, the peak-to-peak total pressure at 
the entrance of the pipe is dependent on the fre- 
quency of rotation of the cam. The data shows a 
droop in the peak-to-peak value of pressure at 
high plenum pressure and high frequencies; these 
trends with plenum pressure and frequency are pre- 
dicted by the computer code. The peak-to-peak to- 
tal pressure at the entrance of the pipe is shown 
to be dependent upon the frequency of rotation of 
the cam. 

Figure 8 displays the peak-to-peak static pre- 
sures versus rotational frequency for various 
plenum pressures for Ute one-per-rev cam. The ex- 
perimental and numeric : curves generally show 

good agreement and the computer code predicts the 
trends. 

Figures 9 and are plots of the ratio of 
pipe end-to- entrance pressure versus cam rotational 
frequency for both cams. The tvo-per-rev cam re- 
sults are included since the two-per-rev cam ex- 
tends the frequency range two-fold. Results are 



Fig. 7 Effect of plenum pressure 



Fig. 8 Peak-to-peak static pressure 

presented for two different plenum pressures for 
both one- and two-per-rev cams. These curves demon- 
strate the presence of a quarter-wavelength reso- 
nance pnenomenon in the pipe. This phenomenon is 
better illustrated in Fig. 9 for the two-per-rev 
cam where a resonance peak occurs at a frequency 
of approximately 44 Hz for the numerical solution 
at the higher pressure and around 55 Hz for the 
corresponding experimental data. The differences 
in che resonant frequencies between the numerical 
and experimental results are likely due to the ex- 
periment having a shorter characteristic pipe length 
because of wave reflection from the end of the noz- 
zle expansion section in the pipe entrance. The 
quasi -one-dimens iona. 1 . model is unable to simulate 
this reflection and instead reflects from the pipe 
entrance itself. It is interesting to ’‘ute that 
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the resonance effect is much more pronounced at 
the higher pressure. For the one-per-rev cam plot 
shown in Fig. 10, no experimental data were obtain- 
ed at the expected resonant frequencies of twice 
the resonant two-per-rev cam frequencies. However, 
the numerical results again show a resonance at 
around 75 Hz for the higher pre: sure and 80 Hz for 
the lower pressure. These resonant frequencies 
are somewhat less than twice the two-per-rev cam 
frequencies. The more pronounced resonant peak for 
the higher pressure in Fig. 10 is consistent with 
the two-per-rev cam results. 



Fig. 9 Resonance phenomenon in pipe; two-per-rev 



Fig, 10 Resonance phenomenon in pipe; one-per-rev 
Effects of Rotational Frequency on Mass Flow Rate 
Figures 11 and 12 show the effect of rotational 


frequency on the computed mean and peak-to-peak 
mass flow rates, respectively, for the one-per-rev 
cam at several different plenum pressures. No ex- 
perimental data are available for comparison with 
the peak-to-peak data. The experimental lines 
shown for the mean mass flow rate represent an ex- 
perimental average since no significant changes in 
average mass flow rate with frequency were obtain- 
ed in the experimental data. Likewise very little 
change was observed in the numerical results. 
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Fig. 11 Mean mass flow rate 



ROTATIONAL FREQUENCY (H*/ 

Fig. 12 Peak-to-peak mass flow rate 


Figure 13 shows the effect of rotational fre- 
quency on the phas? of the first harmonic of the 
mass flow results chown in Figs. 11 and 12. The 
phase shift increases with increasing plenum pre- 
sure (and increasing mass flow rate) and increases 
almost linearly with frequency over most of the 
range examine*!. The linear trend indicates that 
within this range there is a fixed time uelay, r 1 - 
most independent o r frequency, associated with a 
given set of system parameters such as plenum pre- 
sure and blade internal geometry. At the higher 
pressure, capacitance lag is controlling, while at 
the lower rre; sure, the lag should be principally 
sonu, lag, 
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No experimental phase shift data are available 
for comparison because it is very difficult, if not 
impossible, to measure dynamic mass flow containing 
higher harmonics and obtain phase information. 

In summary, the comparison of the numerical 
calculations with experiment, indicate that, even 
with a simple constant discharge coefficient model, 
the computer code satisfactorily predicts the sys- 
tem performance over a wide operating range. It is 
able to predict the trends and magnitudes of the 
total and static pressures as well as the mass flow 
rate for plenum pressures from 0.98 to 13. /5 psig, 
fOi. rotational frequencies from 15 to 120 Hz, and 
for valve area variations from 0 to approximately 
1 in. 


control the rise or decay. Figure 17 shows the pro- 
pagation of the incident pressure wave created by 
sudden valve closing. To smooth these transient 
results, second-order spatial numerical damping has 
been added to the numerical algorithm 1 for computing 
the incremental or '’delta" solution for a given time 
level from the solution for the preceding one. This 
was done by computing the increment from the weighted 
average of the second-order (in time) Beam and 
Warming algorithm and the first-order (in time) al- 
gorithm of Lax\ The weighting factors applied to 
the Beam and Warming method and the Lax method were 
0.98 and 0.02, respe _tively. 



Fig 13 Mass flow rate phase angle 
Theoretical Transient Response 

For pneumatic control system design, in addi- 
tion to a knowledge of the periodic response of the 
system to periodic control inputs, the transient 
dynamics are also of interest. In the present study 
an independent numerical investigation of transient 
dynamics was conducted. This investigation simula- 
ted the pneumodynamic response the model system 
to step inputs in control valve settings. This was 
accomplished by suddenly opening the control valve 
and allowing air to flow into a blade at atmospheric 
pressure, allowing steady flow to be established, 
and then suddenly shutting off the flow, allowing 
the pressure inside the pipe to return to atmospher- 
ic pressure. 

The results of this investigation into the 
system’s transient dynamics are shown in Figs. 14 
to 17. Figure 14 shows the behavior of the mass 
flow rate in response to the two control actions 
for two plenum pressures. The rise or decay time 
lor the mass flow is approximately 0.003 seconds 
and dees not change significantly for the two 
plenum pressures. This suggests that in these 
cases tne characteristic time phenomena are princi- 
pally related to the time it takes for a pressure 
disturbance to propogate the length of the tube, 
since this is also about 0.003 seconds. This ex- 
planation is reinforced by examining Figs. 15 and 
16 which are plots of the pressure distribution zi 
various elapsed times after the sudden valve open- 
ing. Fig. 15 s^ows the propagation cf the incident 
pressure wave created by the valve opening and Fig, 
16 shows the propagation of the wave after its re- 
flection from the end of the tube. The sonic lag 
is apparently controlling the transient response 
to the step input. However, for a system where 
the capacitance lag is much higher or lower than 
the sonic lag, the capacitance 1 g would tend to 



Fig. 14 Theoretical transient mass flow rate 
Conclusions 


Th* results presented indicate that quasi-one- 
dimensional unsteady flow theory can be applied to 
predict, ' ith reasonable accuracy, the pneumodyna- 
mic response of an idealized circulation control 
rotor model to cyclic control valve inputs* 

The results also chow that, for a given set of 
system parameters, the phase lag in the response ' f 
the system to cyclic control input is a fixed time 
delay almost: independent of frequency. Higher har- 
monic content of the mass flow rate output exhibits 
some amplication due to nonlinear effects, but the 
higher harmonic content of the pressure response 
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Fig, 15 Theoretical incident transient pressure 
wave 

does not. 

A resonance phenomenon is predicted by the 
numerical results at approximately the frequency 
for quarter-wave acoustic resonance in a pipe. 

The predicted resonant frequencies are slightly 
different than the experimentally observed resonant 
frequencies. A possible explanation for this is 
the inability of the quasi-one-diraensional theory 
to predict wave reflection from the expansion in 
the entrance section of the experimental model. 

Numerical investigation of the response of 
the flow variables to step inputs Ja the control 
valve area indicate that, for the cases examined, 
the response time is controlled by the character- 
istic time for wave propagation (sonic lag). 
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DYNAMIC RESPONSE CHARACTERISTICS OF A CIRCULATION CONTROL ROTOR MODEL PNEUMATIC SYSTEM 
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Andy Lemnios, Kaman Aerospace : First, an observation. On the bench test that we did on the 

half scale model of the CCR several years ago I believe we observed about a 20 to 30 degree 
pha~e lag between the control input to the control valve to the actual buildup of pressure in 
the plenum chamber. We also observed that the duct of the blade itself charged and discharged 
as a plenum. Basically there was very little phase lag relationship between the root end of the 
blade and the tip end of the blade as far as pressure buildups were concerned. Those pressure 
variations have been documented. I was very pleased to see that in your high freqeuency pneumo- 
dynamic analysis you have included flexibility in there. I was wondering if you have any plans 
to validate the analysis? 

Reader: Yes, Andy, over the years we have done that. We've used the H2 data [that] 1 had to 

a limited extent. The problem there is that all the tests that we've run so far we've always 
found we didn't have all the information we needed. In that rotor system we had total pressure 
at the inboard, [but] the outboard was very questionable. I've had very good results with this 
code with the H-2 experience, that is, if I look at the components I was able to look at the 
weight flow, the one per rev and two per rev component at the inboard station and get good 
correlation if I did a parametric study to identify it. We've looked at the Boeing Vertol data 
which was tested in Boeing's tunnel about two years ago. We've looked at fixed system results 
similar to these; half-scale [tests] that were done at UTRC about 3 or 4 years ago. We have not 
been fortunate enough to get this type of correlation across a large range of parameters, but we 
can always trace it to some unknown that we haven't been able to measure. One of the big prob- 
lems we end up with the flexible slots is identifying what that slot area really is. The 
25-foot rotor system, we ran correlation with that and we came back to the same position. Most 
of the time you end up measuring weight flow which goes into all the blades and you have to 
assume what it is for one blade and then the area variation. I think looking in the box we've 

worked with most of the data that exists and were fairly reasonable based on engineering level. 

Lemnios : Getting back to your comment on the H-2, Ker, the half-scale data I was referring to 

was the one that was done at Honeywell under a subcontract. 

Reader : Yes, I'm aware of that. 

Wayne Johnson : I would like to put Andy's question a little bit differently. You describe what 

you have done in the past and described the capabilities of the analysis that you have. If you 
were making a list of the things that you didn't know about the duct aerodynamics and tha^ 
would really like to know next, where would you start? 

Reader : First of all I think the code gives us a good representation analogy. It is a fairly 

^ood design tool. The [major] problem is we need more instrumentation on the blades; we need 

to know more specifics at different radial stations and that is being done. 

Johnson : So until you do another round of experiments you really don't know . . . 

Reader : That's right. We are using it as a design tool now based on our experiments up to 
now. The next thing we need is more data. I think we will be getting that, probably in the 
next year. 

Johnson : Since I have the microphone I'm going to ask another question. To make a connection 

with some of the things we have heard earlier--the dynamic inflow [models] for all their sim- 
plicity are extraordinarily useful because they are a discrete state model of the dynamics. I 
know, at least I think I've seen in some of the earlier work on circulation control rotors that 
people have developed rather heuristic models of duct aerodynamics, mass-spring type models, 
which I think they Just sort of guessed at Do you think there is a chance that one mig t start 
with the kind of model you are using and then use that as a basis for r more solid foundation 
for a discrete state model for the aerodynamics that then could be a little oit more useful? I 
think that with a two-dimensional Beam-Waiming type solution that y^u are dealing with fairly 
significant computation time, I was lust wondering if there are ways of that sort that you 
might s to make a better, but faster model. 

Ren . There's one way and it was addressed tarlier with Bob Jones' work. If you believe 
thw t'de then you can regress a model that you need to operate with. But the thing of it is 
yoj need to identify all of your parameters up front. There's been some work done on that. 
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That is, you take the code, you can design a system, get it set up, and then you can linearly 
regress out of it a faster working model. One of the problems that we have had historically and 
what has led us to this code is that if you use your analogies to develop [the code] you end up 
with a lot of constraints on it. Now we had that with the SUPERFLOW code, that is, where we 
tailored it for specific configurations so that it’s good for that specific configuration and 
it may run a little bit faster— not a whole lot— a little bit, but you ended up without a good 
design tool. If you go to a real simple model like we had done in the past for stability and 
control then you end up without the higher harmonic capability. So you end up with two codes 
and now you have to start keeping track of two codes. I think you are probably better off using 
a code like this and comt up with a regressed model out of it for your actual application. 

Bob Wood, Hughes Helicopters : I was interested in your model— the traveling wave [going] down 

the blade and then reflected back. I didn't really see the thing coming to equilibrium. In 
other words, is that wave— my question is in several parts, does that wave then come down and 
then travel back? Essentially, what is the amount of time that's required for a step input 
before the flow comes to equilibrium? 

Watkins : In actuality it's going to reflect back and forth any number of times, but the phe- 

nomena tends to be obscured by the flow itself. But essentially it only takes one complete 
reflection for the flow in the duct to pretty much distribute itself evenly over the length. 

Reader : I'd Just like to make a comment on Wayne's question Just a minute ago. I Just had 

another thought. One thing that we f<-jnd, for instance right now, to do this work we've run 
16 cements in it and it has taken the cost up for running the code. However, if you back off 
the number of elements, depending upon what you are using it for, you may not need the fidelity, 
sc if you back off the elements your time goes way down— probably a factor of 4 to maybe 5 
depending upon the number of elements. Kamar. did some of that work earlier with their codes 
and they found that, for instance, for some cases they needed 4 to 6 elements, but for general 
handling quality-type stuff they only needed one. I think that is what Andy was alluding to 
earlier. Depending upon what you are loing with the pneumodynamics and what you need out of 
it— if you only need one per rev cyclic in order to do stability and control you can come up 
with a fairly simple model for that. 
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